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Abstract. We show analytically that, under certain assumptions, the periodic Anderson model and the 
Hubbard model become equivalent within the dynamical mean field theory for quasiparticle weight Z —* 0. 
A scaling relation is derived which is validated numerically using the numerical renormalization group at 
zero temperature and quantum Monte Carlo simulations at finite temperatures. Our results show that the 
/-electrons of the half-filled periodic Anderson model with nearest neighbor hybridization get localized 
at a finite critical interaction strength U c , also at zero temperature. This transition is equivalent to the 
Mott-transition in the Hubbard model. 
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PACS. 71.27.+a Strongly correlated electron systems; heavy fermions - 71.10.Fd Lattice fermion models 
(Hubbard model, etc.) - 71.30.+h Metal-insulator transitions and other electronic transitions 
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The periodic Anderson model (PAM) and the Hubbard 
model (HM) are two of the most fundamental models in 
condensed matter physics. Despite the simplicity of their 
Hamiltonians, the many body nature of these models re- 
sults in a complicated correlated electron problem which 
does not allow for an exact solution except for the one- 
dimensional HM. Also in infinite dimensions where 
both models map onto a single impurity Anderson model 
with different self-consistency conditions, an exact solu- 
tion is only possible numerically. In a recent paper S, 
such a numerical study showed an astonishingly similar 
behavior of the two models and, in particular, that the 
PAM exhibits a transition similar to the Mott transition 
of the HM. 

In the present paper we show analytically that under 
certain assumptions the infinite dimensional PAM and the 
HM become equivalent in the limit of vanishing quasipar- 
ticle weight Z. These assumptions are similar in nature to 
those employed for the self-consistent projective method 
H and the linearized dynamical mean- field theory @. 
Thus, the hitherto numerically found similarity can be un- 
derstood and the critical Coulomb interaction U c for the 
Mott transition of one model can be determined from that 
of the other model. Applying the numerical renormaliza- 
tion group (NRG) at zero temperature (T = 0) and the 
quantum Monte Carlo (QMC) technique at finite T, we in- 
vestigate to what extent the two underlying assumptions 
are fulfilled and show, for the first time, that the zero tem- 



perature PAM with nearest neighbor hybridization has a 
Mott transition with quasiparticle weight Z — > at a fi- 
nite U c , in contrast to the single impurity Anderson model 
for which Z — > at U c = oo. 

The Hamiltonian of the HM reads 



-)<n 



(1) 



Here, fj a and f- a are creation and annihilation operators 
for an electron with spin a on site i or j, respectively, 
n{ a — fl a f i(T , (ij) denotes the sum over nearest neighbors, 
and t the hopping amplitude between them. We use the 
symbol / for the electrons of the HM since the equivalence 
of these and the /-electrons of the PAM will be reported 
in this paper. The PAM consists of a band of conducting 
electrons (d-electrons) and interacting /-electrons. Both 
are coupled via the hybridization Vij : 

H=-tJ2 4 a d ja +J2Vi 3 {d\j ja + h.c.) 



+v y £('4t-b«-b- 



(2) 
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We only consider the particle-hole symmetric case (/i = 
in this form for a symmetric non-interacting density of 
states (DOS). In infinite dimensions or with the number 
of nearest neighbors Z — * oo, a non-trivial scaling of the 
kinetic energy is obtained by t = t* j\[Z. In the following, 
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t* = 1 sets the energy scale. We consider two different 
kinds of hybridizations: (i) a nearest neighbor hybridiza- 
tion Vij — tdf I yf~Z for nearest neighbors i and j which is 
zero otherwise and (ii) an on site hybridization with 
Va — tdf and zero otherwise . 

Within dynamical mean field theory (DMFT) §],§], 
which becomes exact in infinite dimensions, the HM and 
the PAM map onto the same single site problem (which 
depends on the Green function G/, self-energy Ef, T, and 
U) but different self-consistency conditions. For the HM 
this self-consistency condition at frequency lo is given by 



de 



N(e) 



(3) 



where N(e) is the non-interacting DOS. In the case of the 
PAM, the ci-electrons can be integrated out since they en- 
ter only quadratically in the Hamiltonian and the effective 
action. This results in an effective /-electron problem with 
a self-consistency condition that reads 



de 



N(e) 



w-Ef{u)-t%e*/{LJ-e) 



(4) 



for the PAM with nearest neighbor hybridization and 



de 



N(e) 



(5) 



for the PAM with on site hybridization, where N(e) is the 
free d-electron DOS. Note, that the effective one-particle 
potential of the PAM (oc 1/{uj — e)) is frequency depen- 
dent, i.e., retarded, due to the fact that the electrons may 
move from the /-orbitals to the d-band and return at a 
later time. The main difference between nearest neighbor 
[Eq.([|)] and on site hybridization [Eq.(||)] is that the for- 
mer describes metallic f-electrons at U = and within 
a Fermi liquid phase while the latter describes Kondo- 
insulating f-electrons, i.e., a gapped f-electron quasiparti- 
cle peak induced by the hybridization. 

The equivalence of the PAM and the HM (at the Mott 
transition Z — » 0) is shown on the basis of two assump- 
tions: (i) that the metallic phase of the /-electrons may 
be described by Fermi liquid theory at low energies and 
(ii) that the remaining spectral weight of 1 — Z is con- 
tained in two Hubbard bands centered around ±U/2 and, 
in particular, that differences in the internal structure of 
these bands have no influence on the low-energy physics. 
Assumption (ii) is certainly only fulfilled approximately 
and becomes justified if the high energy features are well 
separated from the low-energy features || . For a detailed 
discussion on this assumption see Sec. 2 of || . 

With assumption (i) and Z = (1 — 9Z , //9w| a)=0 ) _1 the 
low-frequency self-consistency condition for the PAM with 
nearest neighbor hybridization is given by 



G f (u)= de 



ZN(e) 



oj - Zt%e*/{u - e) 



(6) 



de- 



ZN(e) 



l+4Zt%- 



u-e/2 + ^Jl + iZtl f e/2 2^1+AZt 



ZN(e) 



l+4Zt%-l 



(7) 

u -e/2-^l+AZi 2 df e/2 2^/1+4^ 
With the partial fraction decomposition above and a vari- 



able transformation y — 2^(1 T yl + 4Zi^)e for the two 
terms of Eqn. ^ one obtains 



dy 



ZN(y) 

w- Zy 



(8) 



which is just the low-frequency self-consistency condition 
of the HM within a Fermi liquid phase. However, with a 
DOS which depends on Z (which is itself a function of U): 



N{y) =N[ 



2Zy 



z 



'l+AZtjj-l) Jl+AZt%-l Jl+AZt% 



2Zy 



MZt%-l 



z 



'l+4Z%+lJ Jl + AZt%. + l JMZt% 



,(9) 



In the limit Z — * 0, N(y) reduces to 



N(y) 



N{y/t%)lt% 



0(Z 2 ). 



(10) 



and becomes, thus, independent of Z. Therefore, at Z — * 
0, the low energy spectral function of the PAM is identi- 
cal to that of the HM with the DOS of Eq. (|To|). With as- 
sumption (ii), i.e., that differences in the internal structure 
of the Hubbard bands have no impact on the low-energy 
physics, the PAM with nearest neighbor hybridization is 
equivalent to a HM which has the DOS of the PAM's d- 
electrons renormalized by the factor l/i;L. Thus, the crit- 
ical Coulomb interaction and temperature for the Mott 
transition of the PAM can be calculated from that of the 
HM via 



and T C PAM =^T C HM 



(11) 



For the PAM with on site hybridization, the same pro- 
ceeding yields for Z — > 0: 



dy 



ZN{-t%/y) t% 
u- Zy y 2 



(12) 



Again, this is equivalent to the low-frequency self-con- 
sistency condition of a HM. While the PAM with near- 
est neighbor hybridization (which has metallic /-electrons 
at small U) maps to a familiar HM, the PAM with on 
site hybridization (which yields a Kondo insulating /- 
electron system) maps to a rather unusual HM. If the free 
d-electron DOS of the PAM N(e) is zero for |e| > D, 
the PAM with on site hybridization is equivalent to a 
HM with a gap of size t%/D in the non-interacting DOS. 
This reflects the Kondo insulating nature of this model. 
Furthermore, the DOS of this unusual HM has not a fi- 
nite bandwidth but tails decaying like 1/y 2 for large en- 
ergies. In particular, the standard deviation of this DOS 
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u/t df 2 

Fig. 1. Quasiparticle weight Z as a function of U/tjj for the 
PAM with nearest neighbor hybridization at a) T — (NRG) 
and b) T/t 2 df = 0.05 (QMC) in comparison to that of the HM. 
The QMC data for the HM are partly from fy}. 

is infinite. For such a DOS, the linearized DMFT || pre- 
dicts a Mott transition at U c — oo and for the HM with 
Lorentzian DOS (a similar DOS without gap) it is known 
that U c = oo Q. Thus, at T = the analytic argument 
suggests U c — oo for the PAM with on site hybridization, 
in agreement with recent NRG results j fTo| . Nevertheless, 
at finite T the transition is very similar to that of the HM 
and the PAM with nearest neighbor hybridization || . This 
is due to the fact that the vanishing of the quasiparticle 
peak at a fixed finite temperature is unaffected by the very 
small energy scale present at T = 0. 

With the approximative nature of assumptionfii) in 
mind we investigate now numerically by NRG pjfl and 
QMC |l3| to what extent the derived scaling relations 
hold. To this end, we calculate Z as a function of 
U/tjf for a Bethe lattice. The results are presented in 
Fig. [y which shows that there is indeed a Mott transi- 
tion Z — ■> at T = in the PAM with nearest neighbor 
hybridization. As in the HM [||[l^,|l6|], the coexistence of 
two solutions is observed at T — (Fig. [j] only contains 
the metallic solution obtained with increasing U). Fi g, p] 
validates, furthermore, that the scaling relation Eqn. (|Ll|) 
holds, at least approximately, even though, the actual val- 
ues of U differ by a factor of 10, both at zero and finite 
temperature (at finite temperature, T/t d f was kept con- 
stant according to Eq.(ll)). 

From Fig. [I] the critical value of U for the Mott tran- 
sition, i.e., the vanishing of Z, is determined fl^| . The 
results are compared to prediction (^lj) in Fig. |^. At small 
tdf, the agreement is very good, while there are notable de- 
viations at larger tdf ■ These deviations can be understood 
from the spectral functions discussed below. 

Fig. ^| shows the disappearance of the central quasipar- 
ticle peak at the Mott transition in the /-electron spec- 
trum of the PAM with nearest neighbor hybridization. At 
the same time, the ci-electron spectral function remains 
finite at the Fermi energy. Thus, despite the Mott transi- 
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Fig. 2. Critical value U c for the Mott transition of the PAM 
with nearest neighbor hybridization at a) T = (circles) and 
b) T/t% = 0.05 (sq uares) compared to t^fU^ (solid line) at 
these temperatures. 
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Fig. 3. Spectral function of the /- and d-electrons for the PAM 
with nearest neighbor hybridization at T = (NRG) and at 
different values of U close to the Mott transition. 

tion of the /-electrons the overall system remains metallic. 
W.r.t. the deviations between U c and prediction (pi]), one 
observes that at t% — 0.2 the quasiparticle resonance is 
well separated from the high-energy Hubbard bands, while 
at tdf = 1 there is additional spectral weight very close to 
the quasiparticle resonance. Thus, assumption (ii) is not 
a good approximation for larger tdf with the consequence 
that the analytic calculation based on assumption (ii) is 
less justified and prediction ( pd| ) less accurate. This ex- 
plains the the ^/-dependence of the deviations in Fig. 2 
which is a priori not clear from the analytic calculation. 

Finally, Fig. 4 shows a comparison between the /- 
spectral functions for the periodic Anderson model with 
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Fig. 4. /-spectral functions of the periodic Anderson model 
with nearest neighbour hybridization compared to the corre- 
sponding Hubbard model according to Eq. (10) at U = 0.25, 



t% = 0.1, and T — Q (NRG). 



nearest neighbour hybridization and the Hubbard model. 
The density of states for the Hubbard model calculation 
is chosen according to Eq. (10). As expected, the results 
show a good agreement in the low-frequency part whereas 
the deviations are more pronounced in the high-frequency 
regime. The agreement for small u> is, however, not per- 
fect, even in the limit Z — * 0, as the derivation of Eq. (10) 
is only approximate. 

In conclusion, we showed analytically that the PAM 
becomes equivalent to the HM at the Mott transition 
Z — > if (i) the low energy physics of both models is 
described by Fermi liquid theory and (ii) the high en- 
ergy features are Hubbard bands which are well separated 
from the low energy quasiparticle peak. This allows to 
calculate the critical interaction U c at which Z — > for 
the PAM from that of the Hubbard model. In particu- 
lar, the PAM with on-site hybridization maps to a Hub- 
bard model with gapped DOS and Lorentzian tails which 
suggests U c — oo, while the PAM with nearest neighbor 
hybridization maps to the Hubbard model with the same 
DOS as the d-electron DOS of the PAM. The latter leads 
to the scaling relation JJ^ AM = t^jU^ M . Numerical cal- 
culations employing NRG and QMC yield that the PAM 
with nearest neighbor hybridization has indeed a Mott 
transition at a finite U c and that the above scaling rela- 
tion is correct for not too large values of t<j/ ■ The similarity 
between both models includes the existence of hysteresis 
for the PAM with nearest neighbor hybridization. 
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